第 43 卷 第 2 期 
2023 年 4 月 


桂林 电子 科技 大 学 学 报 


Journal of Guilin University of Electronic Technology 


月 下! 
Vol. 43, No. 2 
Apr. 2023 


一 种 改进 雅 可 比 算法 的 频 域 临界 采样 图 滤波 器 组 
李威 京 ， 菏 俊 正 


(桂林 电子 科技 大 学 信息 与 通信 学 院 , 广 西 桂林 541004) 


摘要 : 频 域 临 界 采 样 岁 滤波 器 组 需要 对 拉 普 拉 斯 矩阵 进行 特 行 


F 分 解 , 这 导致 了 该 框架 计算 复杂 度 过 高 。 针 对 该 问题 , 采 


用 改进 雅 可 比 算法 近似 求解 该 框架 的 特征 矩阵 ,从 而 降低 计算 复杂 度 。 改 进 的 雅 可 比 算法 将 近似 求解 特征 和 矩阵 的 问题 归 
结 为 一 个 带 约束 的 优化 问题 ,将 拉 普 拉 斯 矩阵 的 近似 误差 作为 目标 函数 ,以 近似 特征 和 矩 阵 的 稀 政 正 交 性 作为 约束 条 件 ,从 


而 求解 出 近似 特征 和 矩阵。 理论 和 仿真 实验 结果 表明 ,近似 特征 矩阵 


j 于 频 域 临 界 采 样 图 滤波 器 组 不 会 改变 其 完全 重 构 条 


件 , 且 与 现 有 的 频 域 临 界 采 样 图 滤波 组 相 比 , 改 进 的 雅 可 比 算法 在 降低 计算 复杂 度 的 同时 保持 了 良好 的 去 噪 性 能 。 
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SApbstract: Critically sampled graph filter banks with spectral domain sampling requires to perform eigendecomposition of the 
aplacian matrix, which leads to high computational complexity. To solve this problem, an improved Jacobi algorithm is 


Proposed to approximate the eigenmatrix of the framework to reduce the computational complexity. In this algorithm, the 


approximate solution of eigenmatrix is formulated into a constrained optimization problem, whose objective function is the 


S approximation error of Laplacian matrix, and the constraint function is the sparse orthogonality of the approximate eigenma- 


+x. Theoretical and simulation experiments show that using the approximate feature matrix in the filter banks will not de- 


sstroy its perfect reconstruction conditions. Compared with the existing critically sampled graph filter banks with spectral do- 


=¥ain sampling, the improved algorithm reduces the computational complexity while maintaining good denoise performance. 
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社交 网 络 、 传 感 器 网 络 及 脑 神 经 网 络 等 一 高 维 
不 规则 数据 可 建 模 为 图 信号 。 与 定义 在 规则 的 时 间 
域 与 空间 域 的 传统 离散 信号 不 同 ,图 信号 通常 定义 在 
不 规则 的 非 欧 几 里 得 域 。 传 统 的 信号 处 理 方法 不 适 
用 于 图 信和 号, 如何 处 理 这 类 高 维 不 规则 数据 已 成 为 一 
个 亟待 解 决 的 问题 。 为 此 ,研究 人 员 提 出 了 图 信号 处 
理 框架 ” ,该 框架 将 传统 信号 处 理 方法 与 图 论 相 结 
合 , 为 处 理 这 类 高 维 不 规则 数据 提供 了 强 有 力 的 工 
有 具 。 其 中 ,由 传统 滤波 器 组 延伸 而 来 的 图 滤波 器 
组 -“ 因 其 稀 玻 特性 和 多 分 辩 分 析 能 力 而 受到 了 广大 
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研究 人 员 的 青睐 。 

目前 ,图 滤波 带 组 的 设计 方法 主要 有 顶点 域 采 样 
图 滤波 器 组 和 频 域 采 样 图 滤波 顺 组 。 文 献 L8] 提 出 了 
完全 重 构 的 两 通道 顶点 域 临 界 采 样 图 滤波 融 组 ,该 滤 
波 融 组 满足 正 交 特性 ,但 不 满足 紧 支 撑 性 。 文 献 L9j] 
提出 了 双 正 交 图 小 波 滤 波 带 组 ,该 滤波 组 既 满 足 完 全 
重 构 特性 ,又 满足 紧 支 撑 性 。 但 上 述 2 种 顶点 域 采样 
滤波 器 组 本 质 上 仅 适 用 于 二 分 图 ,对 于 非 二 分 图 则 需 
要 进行 近似 处 理 。 文 献 L10] 提 出 的 样 条 图 小 波 滤波 
顺 组 满足 完全 重 构 特 性 与 图 不 变性 ,其 对 所 有 拓扑 结 
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构 的 图 信号 都 适用 。 然 而 ,这 几 类 顶点 域 采样 图 滤波 
器 组 还 存在 如 下 局 限 性 :1) 需 要 选择 合适 的 采样 集 ， 
才能 确保 其 完全 重 构 ;2) 其 完 po 
一 ,不 同 的 采样 集会 影响 图 滤波 器 组 的 整体 性 能 
献 [11 提出 oe 
域 采样 图 滤波 器 组 克服 了 顶点 域 采 样 图 滤波 器 的 缺 
点 ,其 完全 重 构 的 采样 集 是 唯一 的 , 且 对 于 任意 拓扑 
结构 的 图 信号 都 满足 完全 重 构 特性 ;但 由 于 其 采样 操 
作 在 频 域 上 进行 ,需要 借助 特征 分 解 来 获取 图 模型 的 
特征 向 量 矩 阵 , 这 导致 了 计算 复杂 度 过 高 。 

针对 频 域 临界 采样 图 滤波 器 组 存在 的 问题 ,可 采 
用 改进 的 雅 可 比 算法 (截断 雅 可 比 算法 和 并 行 截断 雅 
可 比 算法 ) 来 近似 求解 拉 普 拉 斯 矩阵 的 特征 矩阵 。 改 
进 的 雅 可 比 算法 是 一 种 贪 禁 算法 ,其 每 步 迭 代 所 求 得 
的 稀 琉 正 交 矩阵 都 属于 Givens 旋转 矩阵 ,将 每 步 迭 
代 求 得 的 Givens 旋转 和 矩阵 相 乘 ,可 得 到 近似 特征 拢 
陇 (CO 从 而 近似 求 得 图 言 号 的 频 域 表示 。 改 进 雅 可 比 算 
潜在 满足 完全 重 构 的 条 件 下 降低 了 频 域 临 界 采样 图 
庆 焉 器 组 的 计算 复杂 度 。 


图 G 二 (V,E) 是 由 顶点 集 V 二 {0,1,…， 
三 1) 与 边 集 五 组 成 ,其 中 节点 数 N 一 |V | 表示 图 
的 天 小 . 图 G 可 用 权重 和 矩阵 W € R” 表示 ,权重 和 矩 
际 的 元 素 zj 表示 顶点 i 与 顶点 7 之 间 的 关联 程度 ， 
若 w, 二 0， 则 表示 顶点 ; 与 顶点 ) 之 间 没 有 边 相 连 。 
根 骗 权重 矩阵 酌 , 可 定义 度 矩阵 了 . 拉 普 拉 斯 矩阵 工 ， 
度 得 阵 卫 是 对 角 和 矩阵 ,其 对 角 线 元 素 等 于 权重 矩阵 


的 行 和 , 即 a; = wo 而 而 拉 普 拉 斯 矩阵 工 定义 为 


度 矩 阵 与 权重 矩阵 之 差 , 即 工 一 D 一 W。 对 拉 普 拉 斯 
矩阵 工 进行 特征 分 解 , 有 

L=UAU', (1) 
其 中 ,对 角 和 矩阵 A € RY 的 对 角 元 素 由 LL 的 特征 值 
组 成 ,其 对 角 元 素 按 从 小 到 大 的 顺序 排列 , 即 0== A。 所 
Al A "二 AN ， 特征 矩阵 U 二 Lu。 Ul “NI 的 列 
向 量 对 应 于 工 的 特征 向 量 。 

图 信号 可 表示 为 N 维 向 量 f € RY , 即 f==| f。， 
fi1,…,fnaj ， 其 中 第 i 个 元 素 对 应 于 图 上 第 i 个 
节点 的 信号 值 。 与 传统 离散 信号 类 似 ,图 信号 也 有 对 
应 的 频 域 表示 ,图 信号 了 的 图 傅 里 叶 变 换 定 义 为 


fF =U'/, (2) 
道 图 健 里 叶 变 换 定义 为 


f=Uf,， (3) 
式 (2)、(3) 称 为 图 傅 里 叶 变 换 对 。 
图 滤波 可 从 顶点 域 和 频 域 2 个 角度 考虑 。 在 顶 
点 域 ,图 滤波 可 表示 为 拉 普 拉 斯 矩阵 的 多 项 式 , 图 信 
号 的 图 滤波 过 程 可 表示 为 


y=Hf= (Sok? f， (4) 

图 信和 号 在 经 双 过 图 滤波 后 ， 每 个 顶点 的 信号 值 都 是 
其 p-hop 邻 域 的 信号 值 的 线性 组 合 。 在 频 域 ,图 信号 
的 图 滤波 过 程 可 表示 为 


y=Hf=UH AUTfF=UH(A)F, (5) 


其 中 ， 
H(A) =diag(H (40), 五 (CA1)，…， 五 (AND)) 
为 图 滤波 器 的 频率 响应 , 古 (4,) 为 图 频率 的 多 项 式 。 
2. 2 两 通道 频 域 临界 采样 图 滤波 器 组 的 结构 
图 1 为 两 通道 频 域 临界 采样 图 滤波 器 组 的 基本 
结构 5 ,其 中 
H,(A) =diag(H, (40), H, (MA1) 9 H, (An_ 1)) 
为 分 析 滤 波 器 组 中 第 个 滤波 器 的 频率 响应 ; 
Gi(A)=diag(Gi, M0), Gr Ai) ,ne, Gr (AN 1)) 
为 综合 滤波 器 组 中 第 & 个 滤波 器 的 频率 响应 
[I | 和 S il 一 LIN ,J Nz 为 频 域 上 采样 矩 
阵 ; Iys J ws 分 别 为 单位 矩阵 、 反 向 单位 矩阵 ; S ,。 == 
ST 和 S ,二 S10" 为 频 域 下 采样 矩阵 ; 了 为 输入 
信号 ; 了 为 输出 信号 。 
两 通道 频 域 临界 采样 图 滤波 器 组 的 输入 输出 关 
系 为 


;5 d,0 一 


| 二 
‘T=Uo0Go (A)S ,oS oH, CA)Uu 十 (6) 


UG1(A)S ,1S ,HAVUT, 
根据 输入 输出 关系 , 当 子 带 滤 波 器 的 频率 响应 满足 
os 
GDHOAN 二 GODHIOQAN 一 0， 
| 1; 一 0,1,…,N 一 1。 
(7) 
传递 函数 工 一 c2T, 此 时 该 滤波 器 组 满足 完全 重 构 条 
件 55 。 
关于 Gu ) 、 GD) 的 设计 , 现 有 的 设计 方法 
主要 有 正 交 四 和 双 正 交趾 2 种 类 型 。 两 通道 频 域 采 
样 图 滤波 器 组 的 正 交 设 计 取 决 于 昌 , (4;), 而 Ho G4;) 
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分 析 滤波 器 组 综合 滤波 器 组 
UO DS AS (ND Ut 
1 | | Y 六 
f= | 1 >/ 
| H'(4) Si Si » G (4) > Lu 
1 两 通道 频 域 临界 采样 图 滤波 器 组 


的 设计 借鉴 了 经 典 信号 处 理 中 的 设计 方法 “, 即 


Ha) 一 Bi 人 ( 瑟 ) ， (8) 


其 余 滤波 器 互 1 (4;)、Go C4;)、G1(4;) 都 可 由 五 CA) 
求 得 ， 


H(A.;) =H oAXnw ns 
Go(A;) = HN;), 
G1(2;) = HA,), 


Hi:Q;,)++ HiGQN ,1)=c, 
对 字 双 正 交 频 域 采 样 图 滤波 器 组 而 言 , 其 高 通 
网 四 低 通 滤波 器 定义 四, 即 

YT (Ai) =—=Go AN 1) GEC 一 万 CAN 1)， 
(10) 
GPA) .HCL ) 可 通过 频谱 普 分 解 获得 。 该 设计 方法 
次 候 于 经 典 信号 处 理 中 的 双 正 交 小 波 变换 5 。 在 上 
述 暴 件 的 约束 下 , 双 正 交 频 域 采样 图 滤波 器 组 的 完全 
重重 


(9) 
通 滤波 器 


AN 


向 条 件 " 为 
Go DHo A ) Ho NGo MN) 一 c 


2 频 域 临界 采样 图 滤波 器 组 的 快速 实现 方法 


频 域 临界 采样 图 滤波 器 组 需要 计算 图 信号 的 图 
傅 里 叶 变换 ,因此 需 对 拉 普 拉 斯 矩阵 工 进行 特征 分 


(11) 


解 , 求 得 其 特征 和 矩 阵 U, 该 操作 的 计算 复杂 度 为 
O(n*)。 当 图 的 规模 较 小 时 ,可 快速 求 得 其 特征 矩阵 ， 


但 当 图 的 规模 较 大 时 ,计算 特征 矩阵 习 的 代价 会 变 得 
较 大 。 为 降低 求解 特征 矩阵 的 计算 复杂 度 , 可 采用 
改进 的 雅 可 比 算法 来 近似 求解 特征 矩阵 U00 。 

近似 求解 特征 和 矩阵 的 问题 归结 

近似 求解 特征 矩阵 U 的 目的 是 找到 一 个 近似 的 
矩阵 右 , 使 得 以 下 最 优化 问题 取得 最 小 值 ; 
IL—S1S,…S)ASTST :ST 和， 


2.1 


特征 


min 
A,S1 ,Sa 51 


st S$,€S,AED, (12) 


其 中 : S; [a RY ，7 == 1 2 为 稀 琉 的 正 交 和 矩阵 ; 人 
为 对 角 算 阵 ;DD 为 对 角 和 矩阵 集 ;S 为 稀 疏 正 交 算 阵 集 。 
为 了 衡量 近似 特征 矩阵 上 与 特征 矩阵 U 间 的 计 
算 复杂 度 , 给 出 相对 复杂 度 (relative complexity 
gain, 简 称 RCG) 的 定义 5 。 相 对 复杂 度 是 特征 矩阵 
U 的 非 零 元 素 个 数 与 近似 特征 矩阵 辟 中 的 稀 朴 正 交 
因子 S$; 的 非 零 元 素 个 数 之 和 的 比率 , 即 
| |。 
J 


RCG 二 一 一 一 一 一 ， (13) 
> 8; (a 
其 中 ‖ |。 为 算 阵 口 的 0 范 数 。 
2.2 截断 雅 可 比 算 法 
截断 雅 可 比 算法 是 雅 可 比 算法 5 的 改进 , 雅 可 


比 算法 在 近似 对 角 和 矩阵 苞 ; 达到 某 个 精度 后 停止 迭 
代 , 而 截断 雅 可 比 算法 则 预先 设置 了 人 迭代 次 数 。 对 于 
式 (12) 的 最 优化 问题 , 若 将 稀 玻 正 交 和 抑 阵 $ 约束 在 
Givens 旋转 矩阵 集 R。 中 ,并 预先 设 定 稀 玖 正 交 矩阵 
Si 的 个 数 了, 即 可 得 到 截断 雅 可 比 算法 5 。7” 维 向 量 
的 Givens 旋转 固定 其 中 一 2 个 坐标 ,其 余 2 个 坐标 
则 旋转 一 定 的 角度 ,n 维 Givens 旋转 矩阵 可 表示 为 

1 


G0 = oy 有 


1 

(14) 
其 中 :p、g 为 旋转 坐标 ; 9 E [0,2r] 为 旋转 角 ; c = 
显然 ,Givens 旋转 矩阵 由 p 、g 、9 


cos0; 一 sin0 5 。 
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三 个 参数 所 决定 。 截 断 雅 可 比 算法 是 一 种 贪 禁 算法 ， 
其 每 步 迭 代 都 站 在 寻找 一 个 Givens 旋转 矩阵 ,使 得 
代价 函数 下 降 最 快 -六 , 即 

Si min ||S LS | wu， (15) 


其 中 ，| 4 2 为 矩阵 4 的 非 对 角 线 元 素 的 平方 和 。 
式 (15) 的 解 属 于 Givens 旋转 矩阵 ,其 旋转 坐标 


Pq 对 应 于 和 矩阵 L; 绝对 值 最 大 的 元 素 | /加 | 所 在 的 
行 和 列 ,旋转 角 

1 Ls = | | Zz 17 
0=3arctan( | |- (2k 十 1) 了 ,k GE ZL 


截断 雅 可 比 算法 如 下 : 
输入 : 拉 善 拉 斯 矩阵 ,Givens 旋转 次 数 。 
输出 :近似 特征 矩阵 ,近似 特征 值 。 
初始 化 : L) =L，j =1。 
~ 1) 找 出 拉 普 拉 央 外 阵 L; 中 绝对 值 最 大 的 元 素 


2) 确定 元 素 2% 所 在 的 行 坐标 p 和 列 坐 标 g ,并 
元 素 /多 5 ; 


3) 计算 Givens 旋转 角 ， 


7 
01 一 Sarctan( a]， 
2 2 
_ 1 人 Ly) x 
三 Farctan 2 十 1 
"4) 计算 
Err 0 1; 一 


(cos(20, ) 十 二 (多 一 芒 )sin(20)) ， 


= i1€ {1,2); 
《95) 各 err 0 二 err 0， 则 旋转 角 0 三 0 ; 知 
err_01 二 err_0， 则 旋转 角 0 = 二 0，,; 

6) 计算 c = 二 cos9,s 二 sin9 及 Givens 旋转 和 矩阵 
S; 一 Conooi; 


7) 计算 Givens 旋转 后 的 拉 普 拉 斯 矩阵 卫 六 = 


SJL;S;, 若 j 三 本 则 返回 步 又 1) ,继续 六 代 ; 若 7 一 
了 则 终止 迭代 , 且 输 出 近似 特征 和 矩阵 U 二 $1S,…S)， 
近似 特征 值 A ==diagL,jn。 
2.3 并 行 截断 雅 可 比 算法 

式 (13) 给 出 的 相对 复杂 度 增益 只 是 理论 上 的 度 
量 指标 ,在 实际 运算 过 程 中 ,使 用 近似 特征 矩阵 了 运 


算 的 时 间 开 销 与 其 相对 复杂 度 增益 并 不 完全 成 比例 。 
这 是 因为 ,在 编程 环境 中 (如 MATLAB) ,近似 特征 


和 矩阵 品 二 $15,…Sj 与 向 量 相 乘 是 一 个 串 行 过 程 ,而 
特征 矩阵 UU 与 向 量 相 乘 则 是 一 个 并 行 过 程 ”。 因 


此 ,即使 相对 复杂 度 增益 较 大 ,使 用 近似 UU 运算 的 时 
间 也 可 能 相对 较 长 。 

由 截断 雅 可 比 算法 改进 而 来 的 并 行 截断 雅 可 比 
算法 "中 较 好 地 解决 了 上 述 问题 。 在 每 步 迭 代 过 程 
中 ,截断 雅 可 比 算法 只 做 了 1 次 Givens 旋转 ,而 并 行 
截断 雅 可 比 算法 则 可 做 nn/2 次 Givens 旋转 。 因 此 ， 
对 于 需要 进行 次 Givens 旋转 的 近似 过 程 ,并 行 截 
汤 雅 可 比 算法 只 需 选 择 KK = 二 2J /n 个 旋转 因子 ,其 中 
每 个 旋转 因子 S$ 都 是 由 /2 个 不 相交 Givens 旋转 
所 组 成 的 矩阵 ,其 数学 表达 式 可 表示 为 


Ri 
S,=P Pp', (16) 
Ry; 
cos0, — singb, 
R=—(° ) ie {1,2,.,n/2). 
sin 9; cos 0b, 


ee 可 比 算法 类 似 ,并 行 截断 雅 可 比 算法 也 
过 寻找 矩阵 L 中 绝对 值 最 大 的 元 素来 确定 

Ce 旋转 ,不 同 的 是 ,所 选 出 的 z7/2 个 元 素 必须 保 
证 两 两 之 间 不 处 在 同一 行 或 同一 列 。 在 每 步 迭 代 过 
程 中 ,虽然 并 行 截断 雅 可 比 算法 得 到 的 解 不 是 最 优 
解 ,但 其 实际 运行 时 间 要 远 低 于 截断 雅 可 比 算法 1 。 
并 行 截 断 雅 可 比 算法 如 下 : 

输入 ; 拉 普 拉 斯 矩阵 LL,Gives 旋转 次 数 J 每 次 
夫 代 的 Givens 旋转 次 数 1 二 nn /2。 

输出 :近似 特征 矩阵 辟 , 近似 特征 值 A 。 

初始 化 : Li 一 工 兴 一 1。 

1) 取 拉 普 拉 斯 矩阵 工 , 的 上 三 角 元 素 ,并 按 从 大 
到 小 的 顺序 进行 排列 ; 

2) 取 前 n/2 个 最 大 值 元 素 , 且 保证 两 元 素 不 在 
同一 行 或 同一 列 ; 

3) 根据 取出 的 n/2 个 元 素 ,计算 其 对 应 的 Giv- 
ens 旋转 角 ， 


k yk 
| 


0; 二 5 arctan( 元 .| 

4) 计算 c, = 二 cos0,;,s; 二 sin9,, 并 将 其 卦 人 矩阵 
S; 相应 的 位 置 ; 

5) 若 err_g1 达 err_g;, 则 旋转 角 9 二 01; 若 err_ 
qi 过 err_q; , 则 旋转 角 9 二 0,; 

6) 计算 c= 二 cos09,s 二 sin9 及 Givens 旋转 和 矩阵 
Sj 一 Cno; 

7) 计算 Givens 旋转 后 的 拉 普 拉 斯 矩阵 卫 六 一 
SjL;S;，, 若 j 达 靖 则 返回 步 台 1) ,继续 迭代 ;车 j 二 

本 则 终止 迭代 , 且 输 出 近似 特征 矩阵 U 一 $1S,…S)， 

近似 特征 值 A = diagL,n。 


这 
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2.4 计算 复杂 度 分 析 


2.4.1 截断 雅 可 比 算法 

寻找 矩阵 工 ; 中 绝对 值 最 大 的 元 素 是 截断 雅 可 
比 算法 中 计算 复杂 度 最 高 的 操作 ,其 计算 复杂 度 为 
OC IL; 1 ,), 最 坏 情况 下 ,该 操作 的 计算 复杂 度 为 
O(n”)。 Givens 旋转 前 后 的 矩阵 Lj 与 L; 4 仅 有 也、 
d 两 行 与 p .da 两 列 的 元 素 不 相同 ,其 余 位 置 的 元 素 都 
相同 。 若 当前 迭代 步骤 中 所 选 Givens 旋转 坐标 p 、q 
与 前 次 选择 的 不 同 ,计算 复杂 度 就 能 降 到 O(n)。 实 
际 计算 过 程 中 , 绝 大 多 数 情况 是 2 次 选择 的 坐标 不 相 
同 , 因 此 截断 雅 可 比 算法 的 平均 计算 复杂 度 为 Oni 十 
n] ) 。 求 得 近似 特征 和 矩阵 苹 后 ,该 矩阵 与 向 量 或 矩阵 
相 乘 的 计算 复杂 度 为 OCJ )。 


2TT2 并 行 截断 雅 可 比 算法 

OO 对 矩阵 艺 ,元素 进行 排序 是 并 行 截断 雅 可 比 算 
源 目 计算 复杂 度 最 高 的 操作 。 在 每 步 迁 代 中 ,最 多 需 
要 对 n(n 一 1)/2 个 元 素 进行 排序 ,该 操作 的 计算 复 
茶 询 为 OO?log n)。 由 于 并 行 截断 雅 可 比 算法 只 需 
送信 一 O CTJA) 次 ,因此 整体 计算 复杂 度 为 
OR log ”) 。 在 求 得 近似 特征 矩阵 可 后 ,运用 该 年 
和 二 向 量 或 矩阵 相 乘 的 计算 复杂 度 为 OC(J ) 。 并 行 
虑 断 雅 可 比 算法 中 每 个 旋转 因子 S, 包含 了 n/2 个 
GBens 旋转 ,是 一 个 并 行 计算 过 程 。 

泛 为 了 使 上 述 2 种 近似 对 角 化 算法 的 计算 复杂 度 
亿 尘 精确 对 角 化 ， 需要 确保 迭代 次 数 二 O (mn), 通 
常 玩 代 次 数 设 为 J =O(nlog 2) 。 


3 S 仿 真 结果 与 分 析 


3.1 图 信号 去 噪 


将 截断 和 并 行 截断 雅 可 比 算法 与 文献 [8j、[L9j、 
L111 的 算法 在 图 信号 去 品 上 的 性 能 进行 对 比 ,所 有 仿 
真 均 在 相同 环境 下 进行 。 用 GSPBOX 构建 了 
Random sensor、 Swiss roll、Community 三 种 图 ,图 
上 信号 由 图 拉 普 拉 斯 矩阵 前 10% 的 特征 向 量 求 和 组 
成 ,图 2 为 节点 数 为 128 的 3 种 图 信号 。 实 验 中 ,图 
言 号 所 添加 的 噪声 是 均值 为 0、 标 准 差 为 o 的 加 性 高 
斯 噪声 ,图 滤波 器 组 中 的 低 通信 道 的 阔 值 设 为 0, 高 
通信 道 的 阔 值 设 为 3c。 

表 1 一 3 为 不 同 算法 对 含 品 图 信号 去 噪 后 的 信 噪 
比 ,去 噪 结果 是 50 次 随机 实验 结果 的 平均 值 。 通 过 
对 比 实验 结果 ,可 得 如 下 结论 : 


@- Bd Dd 办 
ae 一 K 从 必 和 
委 - | 二 0.6 
a 0— ?a ~“@ lg 0.4 
® a 


/tN 人 所 四 
SS / | i\ -0.6 
ES hd 
bs S 0 ed 0.8 
妇 一 全 
人 一 人 
1.0 
(a) 传感器 网 络 


/ -0.4 


上 


(c) 社交 网 络 


2 3 种 不 同 拓扑 结构 的 图 信号 


1) 在 Givens 旋转 次 数 相同 的 情况 下 ,截断 雅 可 
比 算法 的 去 噪 效果 总 体 上 优 于 并 行 截断 雅 可 比 算 法 。 
这 是 因为 在 每 步 迭 代 过 程 中 ,截断 雅 可 比 算法 获得 的 
解 都 是 局 部 最 优 解 ,而 并 行 截断 雅 可 比 算法 获得 的 解 
并 非 局 部 最 优 解 ,这 导致 了 并 行 截断 雅 可 比 算法 的 收 
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敛 速度 要 慢 于 截断 雅 可 比 算法 。 

2) 随 着 Givens 旋转 次 数 增加 ,截断 雅 可 比 算法 
与 并 行 截断 雅 可 比 算法 的 去 品 效 果 更 接近 于 文献 
L11 ,原因 在 于 , 随 着 和 迭代 次 数 的 增加 ,2 种 算法 求 得 
的 近似 特征 矩阵 不 断 趋 近 于 精准 的 特征 矩阵 。 

3) 在 迭代 次 数 较 低 时 ,截断 和 并 行 截断 雅 可 比 算 


种 信号 混杂 程度 较 高 ,从 采样 图 信号 中 难以 有 效 滤 除 
噪声 信号 。 而 对 于 频 域 采样 图 滤波 器 组 而 言 , 由 于 品 
声 信 息 大 多 分 布 在 高 频 部 分 ,图 信号 与 噪声 信号 从 频 
域 角度 得 到 了 有 效 分 离 。 因 此 , 频 域 采 样 图 滤波 带 
组 ”的 去 噪 效 果 要 优 于 顶点 域 采样 图 滤波 器 组 ， 
而 改进 雅 可 比 算法 是 对 文献 L11j 的 一 种 近似 ,其 去 噪 


法 取得 的 去 噪 效 果 要 优 于 文献 [8-9j]。 在 顶点 域 采 样 ” 效果 本 质 上 仍 优 于 文献 [8-9]。 
图 滤波 器 组 中 ,图 信号 与 噪声 信号 未 得 到 有 效 分 离 ,2 
表 1 sensor 图 信号 在 不 同 算法 下 去 噪 后 的 信 噪 比 dB 
节点 数 为 128 节点 数 为 256 节点 数 为 512 
0 一 1/16 o=1/8 0 一 1/4 0 一 1/16 0 一 1/8 0 一 1/4 0 一 1/16 o=1/8 0 一 1/4 
含 品 信号 14. 07 8. 29 2.07 14.14 8. 21 2. 14 14. 07 8. 05 2. 07 
文献 [8] 12. 31 8.12 4. 62 12:85 8. 68 BLa 12. 90 8.71 5.17 
文献 [9] 12.59 7. 43 3. 22 12. 36 8.08 3. 86 12. 56 7. 55 3. 20 
11] 19. 89 14.13 7. 99 19. 80 13. 90 7. 89 19.77 13. 70 7.79 
Ra G10g 1) 19.04 13. 83 7.97 18. 98 13. 64 178 18. 80 13. 42 7. 69 
(nlog n) 17. 88 13. 42 7.84 18. 35 13. 42 7. 69 17. 94 13. 14 7. 64 
考 酬 (2n 1og ”) 19. 80 14. 06 8.00 19. 56 13. 83 7. 82 19. 57 13. 63 7.73 
R210g n) 19.43 14.02 8.02 19. 45 13.78 7. 81 19. 16 13. 53 7.72 
EN 
©O 表 2 ”roll 图 信号 在 不 同 算法 下 去 噪 后 的 信 噪 比 dB 
AM 二 节点 数 为 128 节点 数 为 256 节点 数 为 512 
© o=1/16 o=1/8 o=1/4 o=1/16 o=1/8 o=1/4 o=1/16 o=1/8 o=1/4 
仿 品 信号 14.33 8.17 2. 26 14. 22 8.17 2. 05 14. 03 8.07 2.07 
广 访 [8] 13.16 8. 63 5.13 12. 46 8. 52 5. 26 12. 66 8. 56 5. 02 
交 称 [9] 12. 00 7. 42 3. 28 12. 64 7. 26 8727 12. 07 7. 68 3. 25 
起 20. 15 14.08 8. 24 20.01 14. 00 7.75 19. 66 13. 67 7.79 
(nlog n) 19.14 13. 83 8.10 19. 19 13.75 7. 65 18. 93 13. 49 7.72 
or n) 18.71 13. 69 8.07 17. 63 12. 98 7. 46 17. 51 12. 98 7. 61 
截断 (2zlog n) 19.93 13.97 8. 20 19. 81 13. 92 7.71 19. 43 13. 47 7.75 
并 行 (2nlog n) 19.70 14.03 8.21 19. 48 13. 84 7.72 19. 00 13. 60 7.73 
表 3 ” community 图 信号 在 不 同 算法 下 去 噪 后 的 信 噪 比 dB 
节点 数 为 128 节点 数 为 256 节点 数 为 512 
1/16 1/8 1/4 1/16 1/8 1/4 1/16 1/8 1/4 
含 品 信 和 号 14. 31 8.07 2.08 14. 20 8. 10 2.20 14.05 8. 06 2. 06 
文献 [8] 13. 50 8.13 3.13 13. 13 7. 60 2. 58 12. 29 7. 33 2. 61 
文献 [9] 13. 86 8. 94 3.76 12. 94 7. 66 2. 76 13. 29 8. 10 3. 04 
文献 [11] 19. 96 13. 82 8.02 19. 84 13. 88 8.01 19. 65 13.77 728 
截断 (nlog n) 17. 26 13.02 7. 82 16. 46 12.73 7. 63 16. 24 12. 52 7. 49 
并 行 (zlog n) 14. 62 10.72 6.91 Ld.13 10. 24 6. 62 13. 07 9. 98 6. 37 
截断 (2zlog n) 18. 83 13.54 7.98 18. 16 13. 37 7. 82 17. 93 13. 20 7. 63 


并 行 (2nlog n) 16. 28 12. 44 7.59 15. 98 


12. 33 {52 15. 45 12. 10 7. 26 
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3.2 图 像 去 噪 


在 相同 噪声 环境 下 对 比 不 同 算法 对 图 像 的 去 噪 
性 能 。 图 像 可 建 模 为 一 张 图 ,文献 [19] 给 出 了 图 像 的 
8 邻 域 图 表示 ,其 像素 对 应 于 图 节点 ,像素 值 则 对 应 
于 图 信号 ,每 个 像素 与 其 上 下 左右 4 个 像素 及 2 条 对 
角 线 上 的 4 个 像素 相连 接 ,其 边 的 权重 值 由 相 邻 节点 
的 像素 值 决定 2 ,如 图 3 所 示 。 

实验 选用 了 如 图 4(a) 所 示 的 硬币 图 像 , 并 以 峰 
值 信 品 比 (PSNR) 和 结构 相似 性 (SSIM) 作 为 图 像 去 
品 的 性 能 评价 指标 。 图 4 为 不 同 算法 去 噪 前 后 的 图 
像 , 表 4、5 分 别 为 图 像 去 噪 前 后 的 峰值 信 噪 比 和 结构 
相似 性 。 从 图 4 和 表 4、5 可 看 出 , 与 文献 L8-9] 相 比 ， 
改进 雅 可 比 算法 的 去 品 效 果 更 好 ; 与 文献 L11] 相 比 ， 
考 品 效果 稍 差 ,但 随 着 Givens 旋转 次 数 的 不 断 增 加 ， 


< 
Be 
~ 
~ 
~ 


a 


未 发 挥 作 用 ,表现 在 去 噪 后 信和 号 信 噪 比 低 于 含 噪 前 ， 
原因 在 于 图 滤波 器 组 滤 掉 的 高 频 部 分 不 仅 含 有 噪声 


改进 雅 可 比 算法 的 去 噪 效果 不 断 趋 近 于 文献 [11] 的 ”信息 ,还 含有 图 像 的 边缘 信息 ,去 噪 过 程 中 原本 的 图 

效果 。 男 外 ,在 低 噪 声 水 平 下 ,部 分 图 滤波 名 组 。 像 信 息 遭 到 破坏 ,而 真正 的 噪声 信息 未 得 到 有 效 处 理 。 
© 
© 
© 
CN 
© 
(qe) 
CN 
已 ) 
CN 
.之 
5 
GS 
ob 
© 

Ce) 文献 [11] (有 截断 (5nlogn) (g) 并 行 (5nlogn) 
图 4 硬币 图 像 在 不 同 算法 下 去 噪 前 后 的 图 像 
表 4 图 像 在 不 同 算法 下 去 噪 前 后 的 PSNR dB 表 5 图 像 在 不 同 算法 下 去 噪 前 后 的 SSIM 

方法 0 一 10 0 一 15 0 一 20 0 一 25 0 一 30 方法 0 一 10 G 一 15 0 一 20 0 一 25 0 一 30 

含 噪 图 像 28.13 24.61 22.12 20.18 18.60 含 噪 图 像 0.9854 0.9677 0.9444 0.9163 0.8859 

文献 [8] 26.73 24.02 22.42 21.35 20.54 文献 [8] 0.9795 0.9617 0.9448 0.9298 0.9166 

文献 [9] 26.76 23.88 22.14 20.97 20.13 文献 [9] 0.9797 0.9607 0.9416 0.9241 0.9099 

文献 [11] 28.31 26.55 25.22 24.01 22.96 文献 [11] 0.9857 0.9787 0.9713 0.9623 0.9528 

截断 (nlog n) 26.75 24.44 23.03 21.96 21.13 截断 logn) ”0.9799 0.9661 0.9537 0.9412 0.9300 

并 行 G1log n) 26.13 23.67 22.12 21.04 20.24 并 行 rlogn) 0.9771 0.9602 0.9443 0.9289 0.9161 

截断 (2nlogn) 27.30 25.26 23.94 22.89 21.99 截断 (2nlogn) 0.9822 0.9717 0.9619 0.9520 0.9418 

并 行 (2nlogn) 26.72 24.52 23.13 22.12 21.31 并 行 (2nlogn) 0.9798 0.9667 0.9546 0.9431 0.9326 
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拉 善 拉 斯 矩阵 的 特征 分 解 是 频 域 临界 采样 图 滤 
波 右 组 中 计算 复杂 度 最 高 的 运算 , 当 图 的 规模 很 大 


时 ,计算 代价 过 大 。 针 对 该 问题 ,使 用 2 种 改进 的 雅 
可 比 算 法 近似 求解 频 域 临界 采样 图 滤波 器 组 中 的 特 
征 和 矩阵 ,从 而 在 保证 完全 重 构 的 前 提 下 达到 降低 计算 
复杂 度 的 目的 。 仿 真实 验 结果 表明 ,在 迭代 次 数 较 低 
时 ,截断 和 并 行 截断 雅 可 比 算法 的 去 噪 性 能 要 优 于 项 
点 域 采 样 图 滤波 需 组 , 且 接 近 于 传统 频 域 临 界 采 样 图 
滤波 顺 组 。 


Ht 
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